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Abstract 



Chaplygin anisotropic matter governing the cosmological evolution in two identical 
universes with an intermediate static spherically symmetric region is considered. The static 



region contains a wormhole allowing one to pass between two horizons. The metric in the 



non-static region represents a kind of an anisotropic Kantowski-Sachs cosmological model 
starting from a horizon instead the initial singularity. A classification of the cosmologies 
^T\ • with a monotonic late stage is presented. It is shown that only one scenario can involve 

a de Sitter regime. The scales in the de Sitter phase allows one to describe the earliest 
accelerated expansion on the classical level. 
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1 Introduction 

One of the fundamental problems of the relativistic cosmology is the problem of the initial 
cosmological singularity. The direct way to solve the problem implies quantization of gravity. 
This line of research develops due to efforts in the string theory, M-theory and loop quantum 
gravity and achievements in high-energy physics also. An alternative view on the origin of 
the Universe supposes to exclude the singularity in its geometrical sense. This approach is 
undoubtedly stimulated by the extensive studying of the modern accelerated expansion of the 
Universe 1 and involves dark energy models such as the cosmological constant, phantom fields, 
exotic equations of state, etc., as well as a modified gravity. The universe evolution without 
the singularity can be realized in some ways. Among them are eternally existing universes 
including emergent universes [2], cyclic evolution [3] and bouncing cosmologies [4]. It is worth 
to mark out the ekpyrotic scenario [5] because of its intriguing confrontation to the inflationary 
paradigm. 

One further feasibility to avoid the initial singularity is universes expanding from a Killing 
horizon [6, 7, 8, 9, 11, 10, 12]. This approach arises from the idea that the universe can 
emerge in an interior of a black hole. More than a half of a century Gliner [13] speculated 
that a de Sitter world can be result of a gravitational collapse. Pathria [14] proposed a more 
definite model of a cosmic expansion beyond an event horizon. The original attention to 
universes inside the black holes was concentrated on the inquiry of geometrical properties and 
cosmological details rather than on exact solutions of the Einstein equations [6] . Models with 
a vacuum source providing a regular black hole core were suggested in [7]. Bronnikov et al. 
[8] used the exact solutions for phantom scalar fields to study and classify regular black holes 
where singularity is replaced by a cosmological expansion and which therefore were called black 
universes. Such objects can be supported by phantom matters and could arise from collapse 
in another ambient universe. In the work [9] the Lemaitre type cosmologies starting from a 
horizon for an anisotropic perfect fluid with the vacuum equation of state were considered for 
different topologies of the spacetime. The authors suggested to call this type of the universe 
origin a null big bang (NBB) since a horizon is a null surface where the spatial volume vanishes. 
A regularity conditions were studied for the NBB cosmologies with neutral and charged matter 
obeying the barotropic equation of state [10], and matter with non- vacuum behavior when the 
equation of state parameter w = p r /p 7^—1 everywhere [11]. 

We construct a similar structure with a Chaplygin wormhole. This kind of wormholes 
can not give an asymptotic flatness and encourages singularities [15] in a completely static 
spacetime. The wormhole surrounded by two horizons naturally avoided these disadvantages. 
The static region connects two cosmological regions and the horizons appear as a NBB. In [12] 
we derived the regularity conditions for the metric and the Chaplygin matter on the horizon 
and obtain some exact solutions corresponding to different cosmological scenarios. 

It was also shown that the considered spacetime configuration is anisotropic in both, static 
and non-static, regions. The observational data are evidence that our Universe is isotropic to a 



1 Note, that many of cosmologies without the singularity was invented long before the discovery of the 
present-day Universe acceleration (see [1] for example) and has been incorporated into dark energy models in 



great extent. This does not forbid, however, an anisotropic behavior at early stage. Our main 
motivation for this paper is to impart more validity to the considered model through conditions 
providing isotropization after an anisotropic phase. We keep it in mind when analysing and 
classifying cosmologies concerned. 

The paper is organized as follows. In the next section we consider the general properties 
of spacetimes with a spherically symmetric wormhole surrounded by two horizons and external 
universes which are supported by matter with the Chaplygin equation of state. The static 
region is described in Sec. 3. The classification of cosmological scenarios is presented in Sec. 4. 
In Sec. 5 we summarize the results and make concluding remarks. 

2 Structure of the spacetime 

We start our consideration with a static spherically symmetric spacetime containing an anisotropic 
perfect fluid in which an energy density and a radial pressure are related by Chaplygin's equa- 
tion of state 

Vr = , (1) 

P 

where a is a positive constant. The metric of the spacetime is taken in the form 

ds 2 = Air) dt 2 tttt - r 2 ( dd 2 + sin 2 9 d<^ 2 ). (2) 

1 — b[r)/r 

Two independent components of the Einstein equations G^ u = 8ttT^ u (we accept c — G — 1) 
give 

V = 8irr 2 p, (3) 

A' b + 8irr 3 p r 

~A ~ r 2 (l - b/r) ' ( ' 

where the prime denotes the derivative with respect to r. Third equation 

Pr = -(Pt-Pr)-(P + Pr)^, (5) 

where pt{r) is the transversal pressure, follows from the energy-momentum conservation law 
V M T MJ/ = providing an integrability of the Einstein equations. 

The metric (2) describes a static spherically symmetric wormhole if the functions A(r) and 
b(r) satisfy the flare-out conditions 

b(r min ) = r min , &'( r min) < 1, (6) 

< A(r min ) < 00, -4'( r min) < 00. (7) 

The radial coordinate r increases from a minimum value r m i n to infinity. Without loss in 
generality we take r min = 1. At this radius the coordinate spheres achieve a minimum area 
and therefore the surface r = 1 is interpreted as a wormhole throat. In this approach r 
simultaneously corresponds to two maps covering the spacetime on either side of the throat. 





Figure 1: The functions A(r) and b(r) for 
the geometry with a wormhole at r = 1 
and simple horizons at r = r . 



Figure 2: The global Carter-Penrose structure 
diagram of the model. The thick lines corre- 
spond to r = oo and the thin lines correspond to 
the horizons at r = r . 



Traversable wormholes also require that A(r) should be positive and b(r) < r elsewhere to 
avoid horizons and singularities. These requirements are ruled out for the considered spacetime 
configuration in favor of the horizon located at the radius r = r > 1. 

The horizon corresponds to simultaneous solutions of the equations b(r ) = r and A(tq) = 
as it depicts in Fig. 1. We consider only simple horizons which imply the following decomposi- 
tion for the metric functions near r 

b(r) = r + b'(r )(r - r ) + o(r - r ), (8) 

A(r) oc (r — r ) + o{r — r ). (9) 

In this case r$ corresponds to the simple horizon separating a static R-region at r < r$ 
from a non-static T-region at r > r . Indeed, the square of the normal to a surface r=const 
in the metric (2) is A = b/r — 1, and A < when 1 < r < tq. In this region the spacetime 
is static (R-region) and r is a spatial coordinate. At r > tq the normal is a spacelike vector 
and A > 0. This is a non-static T-region where r becomes a temporal coordinate. It is verified 
using the Kretschmann scalar that the class of metrics under consideration is regular elsewhere 
[12] including the horizons. The global structure of the spacetime is shown in Fig. 2. It is 
represented as diagrams for de Sitter geometry joined together through the R-regions. 

The regularity conditions on the horizon imply a continuous behavior of the pressures and 
the energy density. When crossing the horizon p and p r interchange their positions in the 
stress-energy tensor and for their continuity the condition p r + p = should be fulfilled on the 
horizon regardless of a matter kind [11]. For the Chaplygin gas it yields 



Note, that the equation of state (1) is precisely the same in the R- and T-regions. It means 
that both regions contain identical matter. 

Compatibility of the Einstein equations near the horizon also leads to the following relation 
for the Chaplygin gas 

b'(r ) = Svrr 2 ^. (11) 

and provides the continuity of the transversal pressure that takes the value p*(ro) = &"(ro)/87rro— 
3a/« on the horizon. 

The equation of state (1) relates the energy density with the radial pressure only. It allows 
to solve the problem predetermining the metric function b(r) satisfying the flare-out and horizon 
conditions. The other functions are found from Eqs. (3)-(5) and produce different cosmological 
scenarios in the T-regions which are presented in Sec. 4. 

3 A wormhole between two horizons 

When r increases from 1 to tq the metric (2) describes a static R-region. To appreciate its 
structure it is convenient to transform the metric (2) to the form 

ds 2 = B(x) dt 2 - d:r 2 -r 2 (:r)(d# 2 + shindy? 2 ), (12) 

where r(x) is an inverse function with respect to the proper length 

x= f dr (13) 

J y/l-b(r)/r 

and B(x) = A(r(x)). 

To give the simple description of this region, consider the function b(r) = r(r — l)(r — r )+r. 
In this case the coordinate sphere radius has the simple form 

r(x) = - ((r + 1) - (r - l)cos:r), (14) 

where x lies in the range [— it, it] and 

r _ r e -(n)+l) 2 (2ro-l)(r2-ro+l)arctan((3r-r„-l)/ v /2+ro-r2)/3r4 v /2+r -r2 

B{X) = B °^ (3r 2 - 2(r - l)r + r + ^W^-Wto-V » ( 15 ) 

where r is taken from (14) and Tq is restricted by 1 < r < 2. 

When x — > ±7T, the function B(x) becomes zero providing the evidence for horizons which 
are spaced at a finite distance apart for a static observer. The function r(x) has a minimum 
at x = implying a wormhole throat, so the R region is a spherically symmetric one with the 
wormhole surrounded by two horizons. 

For an arbitrary function b(r) based on the conditions (6) and (8) the qualitative behaviors 
of the metric functions r(x) and B(x) are similar to (14) and (15), and the typical behavior of 
r(x) and B(x) are depicted in Fig. 3. 
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Figure 3: The metric functions depending on the 
proper length x in the R-region. 

4 Kant owski- Sachs cosmologies 

4.1 General properties 

The metric (2) in the T-region can be rewritten as 

dr] 2 



ds 2 



b(rj)/r] — 1 



A(rj) dp 2 - r] 2 dtt 2 



(16) 



where A (77) = —A (77) is the positive function for 77 > r . We introduce the new coordinates, 77 

instead r, and p instead t, to mark the interchange of roles between the spatial and temporal 

coordinates. The integral 

di] 



t = ± , (17) 

J VKv)/v-i 

defines the physical time r for an observer in the T-region. The Carter-Penrose diagram in 
Fig. 2 shows both the contracting region T_ corresponding to r < and the expanding region 
T + corresponding to r > 0. The metric (16) describes a Kantowski-Sachs (KS) cosmological 
model [16]. The spatial part of the KS universe (i.e. surfaces 77 =const) has the spherical 
symmetry but it is highly anisotropic. A topology of the spatial section is R x § 2 and a 
three-dimensional cylinder expands with two scale factors for the spherical and longitudinal 
directions which depend on the proper time r according to a t (r) = 77 (r) and a r (r) = yl 1//2 (7?(r)) 
correspondingly. 

In all cases the evolution starts with a NBB from a horizon at 77 = r$. There is a coordinate 
singularity which the observer interprets as the beginning of the universe. It starts to expand 
anisotropically with the initial values of the scale factors a r = and a t = Tq. The later 



were obtained in [12]. The present objective is to find a spacetime configuration providing 
isotropization at late times to be of interest in cosmo logical context. 

Since the spherical scale factor at linearly increases with the coordinate time i], asymptotical 
equalizing the scale factors implies that the metric function b(i]) becomes monotonic one from 
some instant so that b(i]) oc i] 13 for rj — > oo. It provides a monotonic behavior of the longitudinal 
scale factor. Before this stage the evolution, in principle, may be arbitrary including, for 
example, cyclic or quasi-cyclic regimes. It is readily seen in Fig. 1 that the rate of growth 
/3 > 1 to avoid singularities in the T-region. All the scenarios can be divided into three types 
according to a value of /3. The scenario of the first kind corresponds to /3 < 3, the second kind 
does /3 > 3 and the third one is realized for /3 = 3. All these types are considered below. 

4.2 Type 1 scenarios: eternal anisotropic expansion 

First we consider the case of /3 < 3. It follows from the equations (1) and (3) that the dominant 
term in the large rj expansion of rhs in the equation (4) is i] 5 ~ 213 . It leads to principally 
exponential growth of the function A(jj) with respect to the coordinate time. The integral (17) 
for large r] gives the power dependence between the physical and coordinate times, r oc t^ 3- ^)/ 2 . 
Due to this asymptotic behavior the universe evolution, having started at some instant in the 
past for the KS observer, lasts for infinite time. All this time the universe anisotropically 
expands. In the remote future the scale factor along the coordinate spheres increases with a 
power of the proper time r while the scale factor in the longitudinal direction grows in the 
mixed manner, power and exponential: 

a r oc t k exp r 4 , a t oc r 2/(3_/3) , (18) 

where the constant k is determined by an explicit form of b(t]). 

Asymptotic behavior of the pressures and the energy density can be estimated from the 
equations (1), (3) and (5). In doing so it should be kept in mind that the energy density 
and the radial pressure interchange their roles when crossing the horizon. In accordance with 
an asymptotic behavior of b{rj) we find different dependencies on the coordinate time which 
transform to the /3-independent behavior for the proper time. From observer's point of view 
the radial pressure decreases in its absolute value as p r oc r~ 2 and the main contribution in the 
energy balance is given by the energy density growing quadratically in the proper time, p oc r 2 , 
and the negative transversal pressure which increases in its absolute value as p t oc r 6 . 

To demonstrate the exact solution let us take the metric function b(i]) in the form (the 
growth rate /3 = 1) 

b(r)) = -(rj- l)(r)-r ) +rj. (19) 

V 

The latter metric function is found as 

A( V ) = A ^^ (r, - q ) q5+q6 (V + go)" exp {|r? 4 + |r? 3 + |r? 2 + g 4 r?} , (20) 




6 


1 







' 


i_ 


4 


1 






p/yfa ^^-^ 


: 


2 


- \ 








: 













- 


2 


- 






Prl^Ta 


- 


4 


- 




^/V* 




: 


6 










- 


1.0 


1.5 


2.0 


2.5 3.0 3.5 


4.0 










r 





Figure 4: The increase of the scale 
factors for the solution specified by 
Eqs. (19) and (20) for r = 1.5. 



Figure 5: The normalized energy density, radial 
and transversal pressures as functions of the ra- 
dial coordinate for the type 1 scenario. The dot- 
ted line separates the R- and T-regions. 



where the constants 



"l V r o + 4r + 4 - rl - 2r 



qo 



r d 



Qi 



2(r - 1) 



d+V 

Q3 = qi (r 2 + 2q 2 Q + q 2 /d + l) , 
q$(2d+l) 



r*d 



?2 = gi(r + l), 



to 



g 4 = gi(r + l)(r 2 + g 2 + i; 
qo(r + 1) 



r 5 d(l - ql) 



2\ ' 



q& 



tqO- - ql) 



(21) 

(22) 
(23) 
(24) 



are positive and determined by the horizon radius r only. Besides, the solution implies q < 1, 
so the function A(if) is zero only at the horizon and positive for rj > tq. 

In this solution the physical time r varies directly as the coordinate one when i] — > oo. At 
late times the asymptotic expansion occurs according to (18) with k, = q 5 and (3 = 1. Figs. 4 
and 5 give a pictorial representation of the scenario considered. 



4.3 Type 2 scenarios: universes with finite lifetime 

The second type of scenarios includes the functions b(i]) with the asymptotical growth rate 
/3 > 3. In this case for large r\ the equation (4) yields A' /A oc — 1/rj from whence it follows 
that A(i]) oc i]^ 1 . It is also easy to see from Eq. (17) that if the function b(i]) — > oo faster than 
rf" for rj — > oo, then the physical time has a final point. These reasons lead to the following 
qualitative pattern of the universe evolution. 

The expansion in the T-region starts with a NBB and the ensuing cosmological evolution 
develops during finite cosmological lifetime from the point of view of the KS observer. Through- 
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Figure 6: The scale factors for the uni- 
verses with finite lifetime. 

peaks at some instant and falls to zero as it is shown in Fig. 6. It easy to see from the asymp- 

— 1/2 

totic dependence of the scale factors that at late stage a r oc a t . Note, that the contraction 
in the longitudinal direction occurs even though both, radial and transversal, pressures are 
negative and the energy density disappears towards the end of the evolution. This is typical 
for the KS model where a spherical acceleration is determined by the radial pressure and the 
spherical Habble parameter H t = at/a t , where the dot denotes the derivative with respect to 
the proper time r, while a longitudinal acceleration also depends on the transversal pressure 
and the energy density according to 



1 tt9 'd t 
+ H 2 t + 2- 



1 2 a r 

l2+ n t ~ — 



0-i 



i7rp r , 



4tt (p-p r + 2p t 



(25) 
(26) 



If the energy density is neglected and the rate of growth of H t is slight so that \p r \ > 2\p t \ + 
Ht/4n, then one observes a decelerated expansion or contraction in the longitudinal direction. 
The metric functions 



A( v ) 



4(3 - 2y/2)r] 4 - 11 + %y/2, 



A Q 



(rj 



V2)t 



(l+^2)(2+v / 2)arctan ( (2»}+2+ v / 2)/v /l4 + 8 v / 2) / \/lA+W^ 



V\ V 



f + (2 + V2)rj + 2 + 3v/2 



(27) 



(2f 



give the example of an exact solution associated with the scenarios of the concerned type (/3 = 4) 
with the horizon radius r$ = 1 + \/2. Fig. 6 depicts the scale factors produced by the metric 
functions (27) and (28). 



4.4 Type 3 scenarios: a quasi de Sitter late stage 

Consider now the case of (5 = 3. At the late stage the relation between the coordinate and 



where the constant 7 depends from an explicit form of b{rj). The scenario for the KS observer 
is following. The universe starts its evolution from the horizon and lasts during for infinite 
time. The asymptotical behavior of the scale factors for the longitudinal and lateral directions 
as functions of the proper time when r — > 00 is 

a r oc exp h r T, a t oc exp h t r, (29) 

where the constants h r and h t are the asymptotic Habble parameters in the longitudinal and 
spherical directions respectively. The constant h t is positive and the constant h r can be positive, 
negative and zero. In accordance with the sign, the initial longitudinal expansion can continue, 
give place to the contraction or settle on a constant value at the remote future. In any case the 
energy density and the pressures tend to constants. This regime is very similar to the de Sitter 
state and, moreover, the anisotropic evolution can isotropize if the metric functions admit the 
solution with h r = ht- 

To study this type of evolution in more details, let us use the metric function 

b{n)=d(ri-l)(n-r )(n-l)+ri, (30) 

where the constants d > and I < 1. Solving the equations (l)-(5) under the conditions 
discussed in the Sec. 2 one obtains the metric function A(rj) in one of two forms 



(rj - r )(rj - l) qi exp {g 3 arctan (q 4 [77 — (1 + r + 0/3])} 
r] (rj 2 - 2(1 + r + l)r]/3d + r + I + r l + l/d) q 



A(n) = A — — ^Pl^^^^'H^ [1/ — {±-r iQ-r lj/omj . . 



° r ( - ) 

A( V ) = A { -^li( v - iy^ v _ q5 yr {v _ % y* (32 ) 

V 

depending on the way that the rational expression in the rhs of (4) is represented as a sum of 
terms with minimal denominators. Under the conditions (7) the constant 



2 + rl - 2r d - r 2 (-2d + ^ + 8 + 4d(r - 1] 



I = ^— r '- (33) 

2d(r -l) l ; 

is expressed by two parameters, r^ and d, as well as all the constants qi, which are not repre- 
sented here in details because of their unhandiness. 

Singularities and additional horizons do not appear in the solutions if I < 1. This condition 
restricts the admitted region for the parameters r$ and d that is shaded in Fig. 7. The dark 
grey corresponds to the solution (31) and the light gray does (32). 

As it is seen from (17) at late stage i] oc exp ydr and hence, h r = 'jy/d/2, h t = yd with 

4 + r 4 _ 4rod + 2 d- 6d 2 - r 2 (-2d + y/r$ + 8 + 4d(r - l) 2 ) 

'' = h -■ (34) 

Various values of r and d bring about different kinds of the quasi de Sitter evolution. These 
are shown in Fig. 7 and the corresponding behavior of the scale factors is depicted in Fig. 8. 

The value 7 = represents the solutions describing freezing of the longitudinal expansion 
with the scale factor a r tending to a constant. The solid line in Fig. 7 corresponds to these 





Figure 7: The parametric diagram represents 
different sorts of the solutions with a quasi de 
Sitter behavior at late stage. The area re- 
stricted by the rectangle is represented in the 
top right panel on an enlarged scale. The detail 
description is in the text. 



Figure 8: The scale factors in type 3 sce- 
narios for r = 1.5. The thick line corre- 
sponds to the spherical scale factor at. The 
thin lines represent the longitudinal scale 
factor a r . 

Line 1: 7 = 0, d ~ 0.41, freezing evolution. 
Line 2: d = 2, 7 < 0, contracting universe. 
Line 3: 7 = 2, d ~ 0.233, isotropization. 
Line 4: d — 0.21, 7 > 2, anisotropic ex- 
pansion. Line 5: d — 0.25, < 7 < 2, 
idem. 



longitudinal direction after the initial expansion. In this region7 < 0. The region under the 
line corresponds to solutions with an infinite expansion. The dotted line in this region is 
specified by 7 = 2 and represents the parameters providing isotropization at late times. The 
parameters under the line generate solutions where the longitudinal expansion is faster then 
the spherical one, and the parameters between the solid and dotted lines give inverse expansion 
rates. Solutions of all these kinds are illustrated in Fig. 8. 



An isotropic regime is realized at the late stage if 7 



2 or, what is the same, 



d 



Art 



4r + 73 - 2r 2 Q - 2r + 1 



(35) 



In this case the universe falls into the de Sitter mode where directional Hubble parameters, 
h r and h t , reduce to the same value H^ = yd and both pressures become equal p r = Pt = 
— 3.£/^ :) /87r. This solution is shown in Fig. 8 by the line 3. It is governed by one independent 
parameter r which dictates the rate of inflation at the last stage. 

Driving parameters can increase in number as well as the scenario can become more com- 
plicated when used different functions biri) with the asymptotic behavior b(n) ~ H^W 3 - The 




Figure 9: The Carter-Penrose diagram 
of the model considered with world lines 
of a motionless KS observer (thick line) 
and photons (dashed lines) coming from 
a contracting universe through a worm- 
hole. 

condition specifying the isotropic regime follows from the equations (4) and (11): 

b'(r ) = SH^rl (36) 

As before, it directly relates the de Sitter phase and the NBB. 

5 Conclusion 

In conclusion we summarize the results obtained and discuss some questions concerned. 

We have studied a spacetime with an anisotropic fluid obeying the Chaplygin equation 
of state. The spacetime consists two T-regions separated by a static spherically symmetric 
R-region with a wormhole providing a traverse between two horizons. In the T-region an 
anisotropic KS cosmology is realized. The peculiar feature of all scenarios is that the evolution 
starts from the horizon appearing as big bang for a KS observer. This singularity is apparent 
since the total spacetime is regular. We have focused our attention on the universes with 
monotonic asymptotic behavior to extract a class of metrics admitting an isotropic late stage. 
We have found three types of such cosmologies. The first one describes an universe with an 
ineradicable anisotropic expansion. The second one represents the contracting universe with 
finite lifetime. The third one in turn includes three sorts of evolution following the longitudinal 
scale factor that can increase, decrease or become constant. The connecting feature for these 
cosmologies is the exponential behavior of the scale factors at late times. The isotropic regime 
occurs in this case only. When b(rj) ~ H^rf the universe is going into the inflation phase with 
the Hubble parameter H^. 

It was shown in [9] that in any universe with a NBB a motionless KS observer can get access 
to prehistoric information from a static region coming with photons and particles which have 
crossed the horizon. It means, as applied to our model, that being in an expanding universe 
T + one can observe not only a static region but receive signals from a contracting universe T_ . 
We illustrate it using the Carter-Penrose diagram in Fig. 9. The thick line gives a motionless 



The proper time r is affine parameter along the geodesic. The dashed lines correspond to 
photon trajectories. The straight line corresponds to the photon having crossed two spatially 
separated horizons and polygonal path does the photon having crossed twice the same horizon 
after reflection in the R-region. 

In closing we estimate a size of the static region for the model with inflation. The asymptotic 
Habble parameter, H^, is related with the horizon radius tq through the equation (35). Turning 
back to the dimensional quantities we obtain H^ = c/\/?>rQ when r$ is close to the wormhole 
throat radius r min , and H^ = c/r when r ^> r min . According to present conceptions the 
Habble parameter for initial inflation takes its value between 10 36 and 10 42 sec -1 . It is easy to 
find that the size of the static region (and the wormhole throat also) falls into the broad range 
from ten to ten billions of the Plank length. This estimation is in agreement with the classical 
approach undertaken in the work. 
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